Quark Matter in a Strong Magnetic Background 

Raoul Gatto and Marco Ruggieri 



Abstract In this chapter, we discuss several aspects of the theory of strong inter- 
actions in presence of a strong magnetic background. In particular, we summarize 
our results on the effect of the magnetic background on chiral symmetry restora- 
tion and deconfinement at finite temperature. Moreover, we compute the magnetic 
susceptibility of the chiral condensate and the quark polarization at zero tempera- 
ture. Our theoretical framework is given by chiral models: the Nambu-Jona-Lasinio 
(NJL), the Polyakov improved NJL (or PNJL) and the Quark-Meson (QM) models. 
We also compare our results with the ones obtained by other groups. 



1 Introduction 

Quantum Chromodynamics (QCD) is the gauge theory of strong interactions. The 
understanding of its vacuum, and how it is modified by a large temperature and/or 
a baryon density, is one of the most intriguing aspects of modern physics. However, 
it is very hard to get a full understanding of its properties, because its most impor- 
tant characteristics, namely chiral symmetry breaking and color confinement, have 
a non-perturbative origin, and the use of perturbative methods is useless. One of 
the best strategies to overcome this problem is offered by Lattice QCD simulations 
at zero chemical potential (see 0] [2] [3] |4] for several examples and see also refer- 
ences therein). At vanishing quark chemical potential, two crossovers take place in 
a broad range of temperatures; one for quark deconfinement, and another one for 
the (approximate) restoration of chiral symmetry. 
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An alternative approach to the physics of strong interactions, which is capable 
to capture some of the non-perturbative properties of the QCD vacuum, is the use 
of models. Among them, we will consider here the Nambu-Jona Lasinio (NJL) 
model 0, see Refs. J6| for reviews. In this gluon-less model, which was inspired 
by the microscopic theory of superconductivity, the QCD interactions are replaced 
by effective interactions, which are built in order to respect the global symmetries 
of QCD. Since gluons are absent in the NJL model, it is not a gauge theory. How- 
ever, it shares the global symmetries of the QCD action; moreover, the parameters 
of the NJL model are fixed to reproduce some phenomenological quantity of the 
QCD vacuum. Therefore, it is the main characteristics of its phase diagram should 
represent, at least qualitatively, those of QCD. 

The other side of the NJL model is that it lacks confinement. The latter, in 
the case of a pure gauge theory, can be described in terms of the center sym- 
metry of the color gauge group and of the Polyakov loop Q, which is an or- 
der parameter for the center symmetry. Motivated by this property, the Polyakov 
extended Nambu-Jona Lasinio model (PNJL model) has been introduced |8] |9), 
in which the concept of statistical confinement replaces that of the true confine- 
ment of QCD, and an effective interaction among the chiral condensate and the 
Polyakov loop is achieved by a covariant coupling of quarks with a background 
temporal gluon field. In the literature there are several studies on the PNJL model, 
see Refs. OS El [12 [IS HI El [13 [13 [IS and references 

therein. 

In this Chapter, we make use of the PNJL model to study the interplay between 
chiral symmetry breaking and deconfinement in a strong magnetic background. 
Moreover, we compute several quantities which are relevant for the phenomenology 
of strong interactions physics in presence of a magnetic background. These topics 
are widely studied in the literature using many theoretical approaches. Lattice stud- 
ies on the response to external magnetic (and chromo-magnetic) fields can be found 
in (26l [27j Ell [29l [30] [66] [671 ■ Previous studies of QCD in magnetic fields, and of 
QCD-like theories as well, can be found in Refs. IEI1E3E3E3III)- Self-consistent 
model calculations of deconfinement pseudo-critical temperature in magnetic field 
have been performed 11361 l73l [38l . 

An important motivation for these kind of studies is phenomenological. In fact, 
strong magnetic fields are produced in non-central heavy ion collisions 13911401 H02l , 
For example, at the energy scale for RHIC it is found eB max w 5nv\\ for collisions at 
the LHC energy scale eB max ps \5m\.. In this case, it has been argued that the non- 
trivial topological structure of thermal QCD, namely the excitation of the strong 
sphalerons 11031 , locally changes the chirality of quarks; this is reflected to event- 
by-event charge separation, a phenomenon which is dubbed Chiral Magnetic Effect 
(CME) |E9l[4T]|42]. The possibility that the CME is observed in heavy ion collision 
experiments has motivated the study of strong interactions in presence of a chirality 
imbalance and a magnetic background, see E2 [104] [1051 [TU6l [T07l [1081 and ref- 
erences therein. An experimental measurement of observables connected to charge 
separation has been reported by the ALICE collaboration in 11091 . It is fair to say 
that realistic simulations of heavy ion collisions show that the magnetic fields have a 
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very short lifetime, and decay before the local equilibrium is reached in the fireball. 
Moreover, the magnetic field is highly inhomogeneous. Furthermore, electric fields 
are produced beside the magnetic fields. Therefore, in order to describe realistically 
hot matter produced in the collisions, one should take care of the aforementioned 
details. However, for simplicity we neglect them, and leave the (much harder) com- 
plete problem to future studies. 

We mainly base the present Chapter on our previous works l60l |6T1 l95ll . We 
firstly discuss chiral symmetry restoration in a strong magnetic background at fi- 
nite temperature, using the PNJL model augmented with the eight-quark interac- 
tion K6l [47l l48l |49l . In this case we also compute the dressed Polyakov loop in a 
magnetic field. The scenario which turns out from our calculations is compatible 
with that of the magnetic catalysis, in which the magnetic field acts as a catalyzer 
for chiral symmetry breaking. Moreover, we discuss on the role of the entanglement 
NJL vertex on the separation between deconfinement and chiral symmetry restora- 
tion in the background field. Finally, we summarize our computation of the magnetic 
susceptibility of the chiral condensate and of quark polarization in a strong magnetic 
background at zero temperature. We base the latter analysis on the Quark-Meson 
(QM) model, which offers the simplest renormalizable extension of the NJL model. 
Throughout this Chapter we consider QCD in the vacuum, that is at zero baryon (as 
well as isospin) chemical potential. Computations at finite chemical potential are 
present in the literature, as we will mention in the main body of the Chapter. 



2 The PNJL model with a magnetic background 

In this Section, we mainly summarize the results obtained in [60, 61). We consider 
quark matter modeled by the following Lagrangian density: 

Jz? = q (if- Dp - mo) q + ga [(qqf + [qij^qf] 

+g%[{w? + {qiy5*q) 2 } 2 , (i) 

which corresponds to the NJL lagrangian with multi-quark interactions [46). The 
covariant derivative embeds the quark coupling to the external magnetic field and 
to the background gluon field as well, as we will see explicitly below. In Eq. (Q~|i, 
q represents a quark field in the fundamental representation of color and flavor (in- 
dices are suppressed for notational simplicity); mo is the bare quark mass, which is 
fixed to reproduce the pion mass in the vacuum, m z = 139 MeV. Our interaction in 
Eq. (QJ consists of a four-quark term, whose coupling g a has inverse mass dimen- 
sion of two, and an eight-quark term, whose coupling constant g% has inverse mass 
dimension of eight. 

We are considering the effect of a strong magnetic background on chiral sym- 
metry restoration as well as deconfinement at finite temperature. We assume the 
magnetic field to be along the positive z— axis; we chose to work in the Landau 
gauge, specified by the vector potential A = (0,Bx, 0). 
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To compute a temperature for the deconfinement crossover, we use the expec- 
tation value of the Polyakov loop, that we denote by L. In order to compute L we 
introduce a static, homogeneous and Euclidean background temporal gluon field, 
Aq = 1A4 = iX a A1, coupled minimally to the quarks via the QCD covariant deriva- 
tive GO. Then 

L=iTr c exp(ij3A a A?) , (2) 

where J3 = 1 /T. In the Polyakov gauge, which is convenient for this study, Aq = 
/A30 + /Ag0 8 ; moreover, we work at zero quark chemical potential, therefore we can 
take L = V from the beginning, which implies A\ = 0. This choice is also motivated 
by the study of IT731 . where it is shown that the paramagnetic contribution of the 
quarks to the thermodynamic potential induces the breaking of the Z3 symmetry, 
favoring the configurations with a real-valued Polyakov loop. 

Besides the Polyakov loop, it is interesting to compute the dressed Polyakov 
loop l43l . In order to define this quantity, we work in a finite Euclidean volume 
with temperature extension j3 = l/T. We take twisted fermion boundary conditions 
along the compact temporal direction, 

q(x,P)=e-*q{x,0), 9G[0,2w], (3) 

while for spatial directions the usual periodic boundary condition is taken. The 
canonical antiperiodic boundary condition for the quantization of fermions at fi- 
nite temperature, is obtained by taking <p = n in the previous equation. The dual 
quark condensate, E„, is defined as 

Z n (m,V) = J o JjjL — {qq) G , (4) 

where n is an integer. The expectation value (-)g denotes the path integral over 
gauge field configurations. An important point is that in the computation of the ex- 
pectation value, the twisted boundary conditions acts only on the fermion determi- 
nant; the gauge fields are taken to be quantized with the canonical periodic boundary 
condition. 

Using a lattice regularization, it has been shown in P31 that Eq. can be ex- 
panded in terms of loops which wind « times along the compact time direction. In 
particular, the case n = 1 is called the dressed Polyakov loop; it corresponds to a 
sum of loops winding just once along the time direction. These correspond to the 
thin Polyakov loop (the loop with shortest length) plus higher order loops, the order 
being proportional to the length of the loop. Each higher order loop is weighed by an 
inverse power of the quark mass. Because of the weight, in the infinite quark mass 
limit only the thin Polyakov loop survives; for this reason, the dressed Polyakov 
loop can be viewed as a mathematical dressing of the thin loop, by virtue of longer 
loops, the latter being more and more important as the quark mass tends to smaller 
values. 

If we denote by z an element of the center of the color gauge group, then 
E„ —> z"E„. It then follows that, under the center of the symmetry group Z3, the 
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dressed Polyakov loop [n = 1) is an order parameter for the center symmetry, with 
the same transformation rule of the thin Polyakov loop. Since the center symmetry 
is spontaneously broken in the deconfinement phase and restored in the confinement 
phase Q (in presence of dynamical quarks, it is only approximately restored), the 
dressed Polyakov loop can be regarded as an order parameter for the confinement- 
deconfinement transition as well. 



2.1 The one-loop quark propagator 

The evaluation of the bulk thermodynamic quantities requires we compute the quan- 
tum effective action of the model. This cannot be done exactly; hence we rely our- 
selves to the one-loop approximation for the partition function, which amounts to 
take the classical contribution plus the fermion determinant. At this level, the effect 
of the strong interactions is to modify the quark mass as follows: 

M = tn Q -2a-4a 3 gs/gl 1 (5) 

where a = g a (qq) ■ As a further simplification, we neglect condensation in the pseu- 
doscalar channel. We notice that the quark mass depends only on the sum of the u 
and d chiral condensates; therefore the mean field quark mass does not depend on 
the flavor. 

To write the one-loop quark propagator in the background of the magnetic field 
we make use of the Leung-Ritus-Wang method 15511 . which allows to expand the 
propagator on the complete and orthonormal set made of the eigenfunctions of a 
charged fermion in a homogeneous and static magnetic field. This is a well known 
procedure, discussed many times in the literature, see for example l56l l57l |59l l58l ; 
therefore it is enough to quote the final result: 

where Ep(x) corresponds to the eigenfunction of a charged fermion in magnetic 
field, and Ep(x) = yo(£p(jc))^yo. In the above equation, 

P=(p + iA4,0,<2yj2k\Q f eB\, Pz ) , (7) 

where k = 0,1,2,. . . labels the k th Landau level, and =S = sign(Qj), with Qf de- 
noting the charge of the flavor /; is a projector in Dirac space which keeps into 
account the degeneracy of the Landau levels; it is given by 



A, = 5«, [£»+5<?,+i + 0»-$3,_i] + (1 - 8 k0 )I , 



(8) 
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where @ > ± are spin projectors and / is the identity matrix in Dirac spinor indices. 
The propagator in Eq. (|6]l has a non-trivial color structure, due to the coupling to the 
background gauge field, see Eq. Q. 

It is useful to write down explicitely the expression of the chiral condensates for 
the flavor / with f = u,d. The chiral condensate is easily computed by taking minus 
the trace of the /-quark propagator. It is easy to show that the following equation 
holds: 

Here, 

tf[L,L,T\p z ,k) = U A -2J^(L,L,T\p z ,k) , (10) 
and jV denotes the statistically confining Fermi distribution function, 

„, - , * l+lLe^f+Le 2 ^ 
tf(L,L,T\p Z7 k) 



l+3uP a f + 3Le 2 P a f +e 3 P a f ' 

(11) 

where 

a} = P : + 2\Q f eB\k + M 2 . (12) 

The first and the second addenda in the r.h.s. of Eq. d27l ) correspond to the vacuum 
and the thermal fluctuations contribution to the chiral condensate, respectively. The 
coefficient J3# = 2 — d^o keeps into account the degeneracy of the Landau levels. 
The vacuum contribution is ultraviolet divergent. In order to regularize it, we adopt 
a smooth regulator U\ , which is more suitable, from the numerical point of view, 
in our model calculation with respect to the hard-cutoff which is used in analogous 
calculations without magnetic field. We chose 

Ua ~ A^ + (pl + 2\Q f eB\kr ' (0) 



2.2 The one-loop thermodynamic potential 

The one-loop thermodynamic potential of quark matter in external fields has been 
discussed in [36, 38], in the case of canonical antiperiodic boundary conditions; fol- 
lowing ETI . it is easy to generalize it to the more general case of twisted boundary 
conditions: 
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go go f=u4 k J-°° 



T I ^Ift/ glog(l+3L e -^- + 3L e - 2 ^- +e - 3 ^- 

f=u,d k 00 

T I r^log(l + 3Z,-^ + 3L e - 2 ^ +e -^ 



(14) 

In the previous equation the arguments of the thermal exponentials are defined as 

& = <0f ( Pz )± i lVj*i, (15) 

with <p defined in Eq. ([3). 

The potential term [L,L, T] in Eq. (fl4l is built by hand in order to reproduce 
the pure gluonic lattice data iflOl . Among several different potential choices l50l we 
adopt the following logarithmic form l9l [T0l . 

% {L T ^ T] =- ^P-LL + b{T) In [1 - 6LL + 4(L 3 + L 3 ) - 3 (LL) 2 ] , (16) 

with three model parameters (one of four is constrained by the Stefan-Boltzmann 
limit), 



a(T) =fl + ai + fl 2 

( 3 

b{T)=b 3 ' 



(17) 



The standard choice of the parameters reads IflOl ; 

a = 3.51, ai = -2.47, 02 = 15.2, £3 = -1.75. (18) 

The parameter 7o in Eq. ( [ToT i sets the deconfinement scale in the pure gauge theory, 
i.e. T c = 270 MeV. 



3 Numerical results 



In this Section, we show our results. The main goal to achieve numerically is the 
solution of the gap equations, 
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Table 1 Parameters of the model for the two choices of the UV-regulator. 

A (MeV) mp (MeV) g (MeV)~ 2 gg (MeV)~ 8 
N = 5 588.657 5^61 5 x KT 6 6x 10~ 22 



This is done by using a globally convergent algorithm with backtrack 13T1 . From the 
very definition of the dressed Polyakov loop, Eq. |@), the twisted boundary condi- 
tion, Eq. (0, must be imposed only in D v . Therefore, we firstly compute the expec- 
tation value of the Polyakov loop and to the chiral condensate, taking <p = n. Then, 
in order to compute the dressed Polyakov loop, we compute the (^-dependent chiral 
condensate using the first of Eq. ( 1191 , keeping the expectation value of the Polyakov 
loop fixed at its value at (p = n IF2TI . 

Following [60, 61 j we report results obtained using the UV-regulator with N = 5. 
As expected, in the other cases no different qualitative results are found; the pa- 
rameter set is specified in Table [T] In the case N = 5, they are obtained by the 
requirements that the vacuum pion mass is m K = 139 MeV, the pion decay constant 
f K = 92.4 MeV and the vacuum chiral condensate (uu) ss (—241 MeV) 3 . In this 
case, the chiral and deconfinement pseudo-critical temperatures at zero magnetic 
field are T* = T P = 175 MeV. 

The main effect of the eight-quark interaction in Eq. (fl} is to lower the pseudo- 
critical temperature of the crossovers. This has been already discussed several times 
in the literature J47l |46], in the context of both the NJL and the PNJL models. 
Therefore, it is not necessary to discuss it further here, while at the same time we 
prefer to stress the results that have not been discussed yet. 

In order to identify the pseudo-critical temperatures, we have define the effective 
susceptibilities as 

dA 



Xa = (m n ) 8 



dT 



A = a,P,Ei. (20) 



Strictly speaking, the quantities defined in the previous equation are not true suscep- 
tibilities. Nevertheless, they allow to represent faithfully the pseudo-critical region, 
that is, the range in temperature in which the various crossovers take place. There- 
fore, for our purposes it is enough to compute these quantities. In Equation (1201) . the 
appropriate power of m K is introduced just for a matter of convenience, in order to 
have a dimensionless quantity; therefore, g = if A = a, E\ , and g = 1 if A = P. 



3.1 Condensates and dressed Polyakov loop 

From now on, we fix JV = 5 unless specified. The results for this case are collected 
in the form of surface plots in Fig. Q] In more detail, in the figure we plot the chi- 
ral condensate S = (cr/2) 1 / 3 , the expectation value of the Polyakov loop, and the 
dressed Polyakov loop E\ , as a function of temperature and magnetic field. 
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Chiral Condensate 




Dressed Polyakov Loop 




Fig. 1 Chiral condensate, Polyakov loop and dressed Polyakov loop as a function of temperature 
and magnetic field, for the case N = 5. From Ref. 1611 . Copyright(2012) by the American Physical 
Society. 
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Fig. 2 Left panel. Chiral condensate S (upper panel), Polyakov loop (middle panel) and L\ (lower 
panel) as a function of temperature, for several values of the applied magnetic field strength, mea- 
sured in units of m\. Right panel. Effective susceptibilities, defined in Eq. )20t . as a function of 
temperature, for several values of eB. Conventions for lines are the same as in the left panel. From 
Ref. 1 61]. Copyright(2012) by the American Physical Society. 



We slice the surface plots in Fig. Q] at fixed value of the magnetic field strength, 
and show the results in Fig. [2] where we plot the chiral condensate (upper panel), 
the Polyakov loop (middle panel) and E\ (lower panel) as a function of temperature, 
for several values of the applied magnetic field strength, measured in units of ra~ . 
In the right panel, we plot fits of the effective susceptibilities in the critical regions, 
as a function of temperature. The fits are obtained from the raw data, using Breit- 
Wigner-like fitting functions. 

The qualitative behavior of the chiral condensate, and of the Polyakov loop as 
well, is similar to that found in a previous study within the PNJL model in the chiral 
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limit [36j. Quantitatively, the main difference with the case of the chiral limit, is that 
in the latter the chiral restoration at large temperature is a true second order phase 
transition (in other model calculations it has been reported that the phase transition 
might become of the first order at very large magnetic field strengths f35l ). On 
the other hand, in the case under investigation, chiral symmetry is always broken 
explicitly because of the bare quark masses; as a consequence, the second order 
phase transition is replaced by a crossover. 

Another aspect is that the Polyakov loop crossover temperature is less sensitive 
to the strength of the magnetic field than the same quantity computed for the chi- 
ral condensate. It is useful, for illustration purpose, to quantify the net shift of the 
pseudo-critical temperatures, for the largest value of magnetic field we have studied, 
eB = I9ni^. In this case, if we take N = 5, then the two crossovers occur simultane- 
ously at eB = 0, at the temperature T* = Tq = 175 MeV; for eB = I9m% we find 
T z = 219 MeV and 7> = 190 MeV. Therefore, the chiral crossover is shifted ap- 
proximately by 25. 1%, to be compared with the more modest shift of the Polyakov 
loop crossover, which is « 8.6%. This aspect will be discussed further in the next 
Section, in which we will comment on the possibility of entanglement between the 
NJL coupling at finite temperature and the Polyakov loop. 

In the lower panels of Figures. Q] and we plot the dressed Polyakov loop as a 
function of temperature, for several values of eB. We have normalized E\ multiply- 
ing the one defined in l43l by the NJL coupling constant. For small values of eB / m^, 
the behavior of E\ as temperature is increased, is qualitatively similar to that at 
eB = 0, which has been discussed within effective models in ll2Tll45l . In particular, 
the dressed Polyakov loop is very small for temperatures below the pseudo-critical 
temperature of the simultaneous crossover. Then, it experiences a crossover in cor- 
respondence of the simultaneous Polyakov loop and chiral condensate crossovers. It 
eventually saturates at very large temperature (for example, in 11211 the saturation oc- 
curs at a temperature of the order of 0.4 GeV, in agreement with the results of B31 ). 
However, we do not push up our numerical calculation to such high temperature, 
because we expect that the effective model in that case is well beyond its range of 
validity. 

As we increase the value of eB, as noticed previously, we observe a tiny splitting 
of the chiral and the Polyakov loop crossovers. Correspondingly, the qualitative be- 
havior of the dressed Polyakov loop changes dramatically: the range of temperature 
in which the E\ crossover takes place is enlarged, if compared to the thin tempera- 
ture interval in which the crossover takes place at the lowest value of eB (compare 
the solid and the dotted lines in Fig. [2] as well as the the lower panel of Fig . [TJ . 

On passing, we discuss briefly the effective susceptibility, dE\/dT, plotted in 
the lower right panel of Fig. [2] since its qualitative behavior is very interesting. 
We observe a double peak structure, which we interpret as the fact that the dressed 
Polyakov loop is capable to feel (and hence, describe) both the crossovers. If we 
were to interpret E\ as the order parameter for deconfinement, and the temperature 
with the largest susceptibility as the crossover pseudo-critical temperature, then we 
obtain almost simultaneous crossover even for very large magnetic field. 
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3.2 Entanglement ofNJL coupling and Polyakov loop 

In l62l it has been shown that the NJL vertex can be deduced under some assump- 
tion from the QCD action; following this derivation a non-local structure of the 
interaction turns out. An analogous conclusion is achieved in l63l . More important 
for our study, the NJL vertex acquires a non-trivial dependence on the phase of the 
Polyakov loop. Therefore, in the model we consider here, it is important to keep into 
account this dependence. Here we follow the phenomenological ansatz introduced 
in ]53], that is 

G = g a [l-a 1 LL-a 2 {L 3 +L 3 )] , (21) 

and we take L = L. Moreover, we mainly discuss here the case without 8-quark inter- 
action. The model with coupling constant specified in Eq. (fJTJ is named Entangled- 
Polyakov improved-NJL model (EPNJL in the following) [53], since the vertex de- 
scribes an entanglement between Polyakov loop and the interaction responsible for 
chiral symmetry breaking. 

The functional form in the above equation is constrained by C and extended Z3 
symmetry. We refer to 1331 for a more detailed discussion. The numerical values of 
OL\ and a-i have been fixed in 1531 by a best fit of the available Lattice data at zero 
and imaginary chemical potential of Ref. l64l . which have been confirmed recently 
in l65l . In particular, the fitted data are the critical temperature at zero chemical 
potential, and the dependence of the Roberge-Weiss endpoint on the bare quark 
mass. 

The values tt\ = 0,2 = a = 0.2 ± 0.05 have been obtained in lf53l using a hard 
cutoff regularization scheme. We will focus mainly on the case a = 0.2 as in l53l . 
In l60l we have verified that in the regularization scheme with the smooth cutoff, the 
results are in quantitative agreement with those of [53 1. There, a detailed discussion 
of the role of a can be found as well (we will skip this discussion in this Chapter). 

We plot in Fig [3] the chiral condensates of u and d quarks as a function of tem- 
perature, at eB = I5m\ and eB = 30m^. In the lower panel of the figure, we plot the 
expectation value of the Polyakov loop as a function of temperature. The conden- 
sates are measured in units of their value at zero magnetic field and zero tempera- 
ture, namely ojj = {uu) = (dd) = (—253 MeV) 3 . They are computed by a two-step 
procedure: firstly we find the values of a and L that minimize the thermodynamic 
potential; then, we make use of Eq. d27l i to compute the expectation values of uu 
and dd in magnetic field. If we measure the strength of the crossover by the value of 
the peak of \da/dT\, it is obvious from the Figure that the chiral crossover becomes 
stronger and stronger as the strength of the magnetic field is increased, in agreement 
with Eg). 

The results in Fig [3] show that, identifying the deconfinement crossover with the 
temperature Ti at which dL/dT is maximum, and the chiral crossover with the tem- 
perature T x at which \do/dT\ is maximum, the two temperatures are very close also 
in a strong magnetic field. From the model point of view, it is easy to understand 
why deconfinement and chiral symmetry restoration are entangled also in strong 
magnetic field. As a matter of fact, using the data shown in Fig [3] it is possible to 
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Fig. 3 Upper panel: Chiral condensates of u and d quarks as functions of temperatures in the 
pseudo-critical region, at eB = 15m^ and eB = 30m^. Condensates are measured in units of their 
value at zero magnetic field and zero temperature, namely o~o = (—253 MeV) 3 . Lower panel: 
Polyakov loop expectation value as a function of temperature, at eB = I5m^ and eB = 30m|. Data 
correspond to a = 0.2. From Ref. 1611 . Copyright(2012) by American Physical Society. 



compute the NJL coupling constant in the pseudo-critical region, which turns out 
to decrease of the 15% as a consequence of the deconfinement crossover. There- 
fore, the strength of the interaction responsible for the spontaneous chiral symmetry 
breaking is strongly affected by the deconfinement, with the obvious consequence 
that the numerical value of the chiral condensate drops down and the chiral crossover 
takes place. We have verified that the picture remains qualitatively and quantitatively 
unchanged if we perform a calculation at eB = 30m^. In this case, we find Tl = 224 
MeV and T x = 225 MeV. 

This result can be compared with the previous calculations ll36l . described also in 
the previous Section, in which the Polyakov loop dependence of the NJL coupling 
constant was not included. In |f36l we worked in the chiral limit and we observed that 
the Polyakov loop crossover in the PNJL model is almost insensitive to the magnetic 
field; on the other hand, the chiral phase transition temperature was found to be very 
sensitive to the strength of the applied magnetic field, in agreement with the well 
known magnetic catalysis scenario |3ll. This model prediction has been confirmed 
within the Polyakov extended quark-meson model in J73], when the contribution 
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from the vacuum fermion fluctuations to the energy density is kept into account Q; 
we then obtained a similar result in AMI , in which we turned from the chiral to the 
physical limit at which m K = 139 MeV, and introduced the 8-quark term as well 
(PNJLg model, according to the nomenclature of ll53l ). The comparison with the 
results of the PNJLg model of [[6T1 is interesting because the model considered there, 
was tuned in order to reproduce the Lattice data at zero and imaginary chemical 
potential [20|, like the model we use in this study. Therefore, they share the property 
of describing the QCD thermodynamics at zero and imaginary chemical potential; 
it is therefore instructive to compare their predictions at finite eB. 

For concreteness, in ||6D we found 7> = 185 MeV and T x = 208 MeV at 
eB = 19m^, corresponding to a split of w 12%. On the other hand, in the present 
calculation we measure a split of ss 1.5% at the largest value of eB considered. 
Therefore, the results of the two models are in slight quantitative disagreement; this 
disagreement is then reflected in a slightly different phase diagram. We will draw 
the phase diagram of the two models in a next Section; however, since now it is easy 
to understand what the main difference consists in: the PNJLs model predicts some 
window in the eB — T plane in which chiral symmetry is still broken by a chiral 
condensate, but deconfinement already took place. In the case of the EPNJL model, 
this window is shrunk to a very small one, because of the entanglement of the two 
crossovers at finite eB. On the other hand, it is worth to stress that the two mod- 
els share an important qualitative feature: both chiral restoration and deconfinement 
temperatures are enhanced by a strong magnetic field; the latter conclusion is in 
qualitative agreement with the Lattice data of D'Elia et al. J26), but in disagreement 
with the more recent data of the Wuppertal-Budapest group ll66ll67l . We will come 
back to a comparison with Lattice data, as well as with other computations, in the 
next Section. 



4 Phase diagram in the eB — T plane 

In Fig. [4] we collect our data on the pseudo-critical temperatures for deconfinement 
and chiral symmetry restoration, in the form of a phase diagram in the eB — T plane. 
In the upper panel we show the results obtained within the EPNJL model; in the 
lower panel, we plot the results of the PNJLs model, that are obtained using the fit- 
ting functions computed in l6D . In the figure, the magnetic field is measured in units 
of m\\ temperature is measured in units of the deconfinement pseudo-critical tem- 
perature at zero magnetic field, namely 7g = o = 185.5 MeV for the EPNJL model, 
and 7g = o = 175 MeV for the PNJLs model. For any value of eB, we identify the 
pseudo-critical temperature with the peak of the effective susceptibility. 



1 If the vacuum corrections are neglected, the deconfinement and chiral crossovers are found to 
be coincident even in very strong magnetic fields 1731 . but the critical temperature decreases as a 
function of eB; this scenario is very interesting theoretically, but it seems excluded from the recent 
Lattice simulations 1261 . 
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Fig. 4 Upper panel: Phase diagram in the eB — T plane for the EPNJL model. Temperatures 
on the vertical axis are measured in units of the pseudo-critical temperature for deconfinement at 
eB — 0, namely T c = 185.5 MeV. Lower panel: Phase diagram in the eB — T plane for the PNJLg 
model. Temperatures on the vertical axis are measured in units of the pseudo-critical temperature 
for deconfinement at eB = 0, namely T c = 175 MeV. In both the phase diagrams, T x , Tl correspond 
to the chiral and deconfinement pseudo-critical temperatures, respectively. The grey shaded region 
denotes the portion of phase diagram in which hot quark matter is deconfined and chiral symmetry 
is still broken spontaneously. From Ref. 1611 . Copyright(20 1 2) by American Physical Society. 



It should be kept in mind, however, that the definition of a pseudo-critical tem- 
perature in this case is not unique, because of the crossover nature of the phenomena 
that we describe. Other satisfactory definitions include the temperature at which the 
order parameter reaches one half of its asymptotic value (which corresponds to the 
T —> limit for the chiral condensate, and to the T —> +°° for the Polyakov loop), 
and the position of the peak in the true susceptibilities. The expectation is that the 
critical temperatures computed in these different ways differ from each other only of 
few percent. This can be confirmed concretely using the data in Fig.[3]at eB = 30m^. 
Using the peak of the effective susceptibility we find T z = 225 MeV and Tl = 224 
MeV; on the other hand, using the half-value criterion, we find T% = 227 MeV and 
Ti = 222 MeV, in very good agreement with the previous estimate. Therefore, the 
qualitative picture that we derive within our simple calculational scheme, namely 
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the entanglement of the two crossovers in a strong magnetic field, should not be 
affected by using different definitions of the critical temperatures. 

Firstly we focus on the phase diagram of the EPNJL model. In the upper panel of 
Fig. [4] the dashed and dot-dashed lines correspond to the deconfinement and chiral 
symmetry restoration pseudo-critical temperatures, respectively. As a consequence 
of the entanglement, the two crossovers stay closed also in very strong magnetic 
field, as we have already discussed in the previous Section. The grey region in the 
Figure denotes a phase in which quark matter is (statistically) deconfined, but chiral 
symmetry is still broken. According to [68 [69), we can call this phase Constituent 
Quark Phase (CQP). 

On the lower panel of Fig. [4] we have drawn the phase diagram for the PNJLg 
model based on Ref. (6l | and discussed in the previous Section. The most astonish- 
ing feature of the phase diagram of the PNJLg model is the entity of the split among 
the deconfinement and the chiral restoration crossover. The difference with the re- 
sult of the EPNJL model is that in the former, the entanglement with the Polyakov 
loop is neglected in the NJL coupling constant. As we have already mentioned in 
the previous Section, the maximum amount of split that we find within the EPNJL 
model, at the largest value of magnetic field considered here, is of the order of 2%; 
this number has to be compared with the split at eB = 2Qm\ in the PNJLg model, 
namely w 12%. The larger split causes a considerable portion of the phase diagram 
to be occupied by the CQP. 



4.1 Comparison with other computations 

In this Section we summarize the main results obtained in the literature, comparing 
them with the scenario depicted in our works. Before going ahead, it is useful to 
summarize the two main results obtained within our one-loop computations: 

• The critical temperature for chiral symmetry restoration is increased by an exter- 
nal magnetic field; 

• The split between deconfiment and chiral symmetry restoration temperatures in 
a strong magnetic background can be reduced if the entanglement vertex is con- 
sidered. 

The first conclusion is in agreement with most of the computations: calculations 
based on the quark-meson-model with and without quantum fluctuations ll72l I7T1 
[73] [75] [78], on chiral perturbation theory at finite temperature l70l . on the PNJL 
model (36] [75), on the holographic correspondence 1791 [80118111771 . 

Besides models predictions, lattice computations of the critical temperatures in a 
magnetic background have been performed 126116611571 . The computations of |[251 
have been performed with a pion mass of the order of 400 MeV, hence a little bit far 
from the physical limit. In this case, both the chiral symmetry restoration and the 
deconfimenent temperatures are measured, and they are found to increase slightly 
with the magnetic field strength; moreover, the two transitions seem to be entangled 
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even at the largest value of the magnetic field considered in the study. On the other 
hand, in ||66l [67) quark masses are chosen such that the lattice pion has its physi- 
cal mass; in this case a non-trivial dependence of the critical temperature for chiral 
symmetry restoration on the magnetic field strength and the quark mass is found. In 
more detail, for the up and down quark condensates at the physical limit, the crit- 
ical temperature decreases with the magnetic field strength. On the other hand, the 
critical temperature for the strange quark condensate is increased by the magnetic 
field, in agreement with the magnetic catalysis scenario. A measurement of the de- 
confinement temperature has not been performed. The results of 1 66 67 1 are quite 
surprising since they reveal an unexpected role of the quark mass on the curvature 
of the critical line in the eB — T plane. 

Computations within the MIT bag model 174] do not have direct access to the chi- 
ral symmetry restoration, but to the deconfinement temperature. Within this model 
it is found that the critical temperature for deconfinement is a decreasing function of 
the external magnetic field strength. This conclusion is in agreement with a previous 
computation [35). Furthermore, a decreasing temperature is found also within the 
quark-meson model if the fermion vacuum contribution is neglected (73) . If decon- 
finement and chiral symmetry restoration are entangled at finite magnetic field, then 
the MIT model based study would give some hint on the mechanism which makes 
the temperature for chiral symmetry restoration in a magnetic background lower 
than that at zero field. However, if this is the case, then the role of the quark mass on 
the dependence of the critical temperature on the magnetic field should be transpar- 
ent. In our opinion, more study is needed to understand the puzzling behavior of T c 
as a function of the magnetic field found on the Lattice: beside model computations, 
independent Lattice simulations should be performed, in order to confirm the results 

of GSIEa. 



5 Polarization of the quark condensate 

It has been realized that external fields can induce QCD condensates that are absent 
otherwise [82]. Here we focus on the magnetic moment, (fE^ v f) where / denotes 
the fermion field of the flavor /-th, and L^ v = -i(y' 1 Y v - y v y^)/2. At small fields 
one can write, according to 182) . 

(f^ v f)=X(ff)Qf\eB\, (22) 

and x is a constant independent of flavor, which is dubbed magnetic susceptibility 
of the quark condensate. In ll82l it is proved that the role of the condensate (1221 
to QCD sum rules in external fields is significant, and it cannot be ignored. The 
quantity % has been computed by means of special sum rules [82, 83 84 85, 86 1, 
OPE combined with Pion Dominance ll87) . holography (88, 89], instanton vacuum 
model [ 90 1 , analytically from the zero mode of the Dirac operator in the background 
of a SI/ (2) instanton 1191) , and on the Lattice in two color quenched simulations at 
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zero and finite temperature |92l . It has also been suggested that in the photoproduc- 
tion of lepton pairs, the interference of the Drell-Yan amplitude with the amplitude 
of a process where the photon couples to quarks through its chiral-odd distribu- 
tion amplitude, which is normalized to the magnetic susceptibility of the QCD vac- 
uum, is possible ||93l . This interference allows in principle to access the chiral odd 
transversity parton distribution in the proton. Therefore, this quantity is interesting 
both theoretically and phenomenologically. The several estimates, that we briefly 
review in Section III, lead to the numerical value of % as follows: 

Xiff) =40-70 MeV. (23) 

A second quantity, which embeds non-linear effects at large fields, is the polar- 
ization, fif, defined as 

E f = {fE 12 f), (24) 

which has been computed on the Lattice in l92l for a wide range of magnetic fields, 
in the framework of two-color QCD with quenched fermions. At small fields [if = 
\xQ/eB\ naturally; at large fields, non-linear effects dominate and an interesting 
saturation of /X/ to the asymptotic value jj.^ — 1 is measured. According to ||92| the 
behavior of the polarization as a function of eB in the whole range examined, can be 
described by a simple inverse tangent function. Besides, magnetization of the QCD 
vacuum has been computed in the strong field limit in |94l using perturbative QCD, 
where it is found it grows as B\ogB. 

In l95l we compute the magnetic susceptibility of the quark condensate by means 
of the NJL and the QM models. This study is interesting because in the chiral mod- 
els, it is possible to compute self-consistently the numerical values of the conden- 
sates as a function of eB, once the parameters are fixed to reproduce some charac- 
teristic of the QCD vacuum. We firstly perform a numerical study of the problem, 
which is then complemented by some analytic estimate of the same quantity within 
the renormalized QM model. Moreover, we compute the polarization of quarks at 
small as well as large fields, both numerically and analytically. In agreement with 
the Lattice results l92l . we also measure a saturation of [if to one at large fields, in 
the case of the effective models. Our results push towards the interpretation of the 
saturation as a non-artifact of the Lattice. On the contrary, we can offer a simple 
physical understanding of this behavior, in terms of lowest Landau level dominance 
of the chiral condensate. As a matter of fact, using the simple equations of the mod- 
els for the chiral condensate and for the magnetic moment, we can show that at large 
magnetic field [if has to saturate to one, because in this limit the higher Landau lev- 
els are expelled from the chiral condensate; as a consequence, the ratio of the two 
approaches one asymptotically. 

We also obtain a saturation of the polarization within the renormalized QM 
model. There are some differences, however, in comparison with the results of the 
non-renormalized models. In the former case, the asymptotic value of [if is charge- 
dependent; moreover, the interpretation of the saturation as a lowest Landau level 
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(LLL) dominance is not straightforward, because the renormalized contribution of 
the higher Landau levels is important in the chiral condensate, even in the limit of 
very strong fields. It is possible that the results obtained within the renormalized 
model are a little bit far from true QCD. As a matter of fact, in the renormalized 
model we assume that the quark self-energy is independent on momentum; thus, 
when we take the limit of infinite quark momentum in the gap equation, and absorb 
the ultraviolet divergences by means of counterterms and renormalization condi- 
tions, we implicitly assume that that quark mass at large momenta is equal to its 
value at zero momentum. We know that this is not true, see for example l97l |%1 : 
even in the renormalized theory, the quark self-energy naturally cuts off the large 
momenta, leading to LLL dominance in the traces of quark propagator which are 
relevant for our study. Nevertheless, it is worth to study this problem within the 
renormalized QM model in its simplest version, because it helps to understand the 
structure of this theory under the influence of a strong magnetic field. 

In our calculations we neglect, for simplicity, the possible condensation of 
p— mesons at strong fields li98l 199] ■ Vector meson dominance [|98l and the Sakai- 
Sugimoto model II 1001 suggest for the condensation a critical value of eB c m mj, m 
0.57 GeV 2 , where nip is the p — meson mass in the vacuum. Beside these, a NJL- 
based calculation within the lowest Landau level (LLL) approximation ll99ll predicts 
p— meson condensation at strong fields as well, even if in the latter case it is hard 
to estimate exactly eB c , mainly because of the uncertainty of the parameters of the 
model. It would certainly be interesting to address this problem within our calcula- 
tions, in which not only the LLL but also the higher Landau levels are considered, 
and in which the spontaneous breaking of chiral symmetry is kept into account self- 
consistently. However, this would complicate significantly the calculational setup. 
Therefore, for simplicity we leave this issue to a future project. 

In 1 95 1 the computation of the polarization and of the magnetic susceptibility 
has been performed both within the NJL and the Quark-Meson (QM) models; the 
qualitative picture does not depend on the model considered. Moreover, within the 
QM model an analytical computation of the aforementioned quantities within the 
renormalized quantum effective potential is feasible (the NJL model, at least with 
a contact interaction as considered in ||95l , is not renormalizable). Therefore in this 
review chapter we limit ourselves to summarize the results obtained within the QM 
model, both numerical and analytical, deferring to the original reference for further 
details. 



5.1 Non-renormalized Quark-Meson model results 

In the QM model, a meson sector described by the linear sigma model lagrangian, 
is coupled to quarks via a Yukawa-type interaction. The model is renormalizable 
in D = 3 + 1 dimensions. However, since we adopt the point of view of it as an 
effective description of QCD, it is not necessary to use the renormalized version of 
the model itself. On the contrary, it is enough to fix an ultraviolet scale to cutoff the 
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divergent expectation values; the UV scale is then chosen phenomenologically, by 
requiring that the numerical value of the chiral condensate in the vacuum obtained 
within the model, is consistent with the results obtained from the sum rules H 1 1 1 . 
This is a rough approximation of the QCD effective quark mass, which smoothly 
decays at large momenta |96ll97l . In Section IV we will use a renormalized version 
of the model, to derive semi-analytically some results in the two regimes of weak 
and strong fields. 

The lagrangian density of the model is given by 



In the above equation, q corresponds to a quark field in the fundamental repre- 
sentation of color group SU(3) and flavor group SU (2); the covariant derivative, 
D^i = <9|i — QfeAp, describes the coupling to the background magnetic field, where 
Qf denotes the charge of the flavor /. Besides, a, % correspond to the scalar singlet 
and the pseudo-scalar iso-triplet fields, respectively. The potential U describes tree- 
level interactions among the meson fields. In this article, we take its analytic form 
as 



where the first addendum is chiral invariant; the second one describes a soft explicit 
breaking of chiral symmetry, and it is thus responsible for the non-zero value of 
the pion mass. For h = 0, the interaction terms of the model are invariant under 
SU(2) V ®SU(2) A <g) C/(l)y. This group is broken explicitly to U(l)l ® U(l)\ <g> 
U(l)v if the magnetic field is coupled to the quarks, because of the different electric 
charge of u and d quarks. Here, the superscript 3 in the V and A groups denotes the 
transformations generated by T3, T375 respectively. Therefore, the chiral group in 
presence of a magnetic field is U(l)y <g> U(l)\. This group is then explicitly broken 
by the h- term to t/(l)y. 

The formalism which is used to compute the magnetic susceptibility and the 
polarization of the quark condensate is similar to the one described in the previous 
Sections; therefore it is not necessary to give the details here. It is enough to write 
down the expresisons for the chiral condensate at zero temperature, 



where the divergent integral on the r.h.s. of the above equation has to be understood 
regularized as in ( fT3l ), and for the magnetic moment for the flavor /, 



if = q [iDpf- - g(a + iyst- &)] q 




(25) 




(26) 




(27) 




(28) 



We take B = (0,0, B); in this case, only E 12 =•£/ is non-vanishing. Using the prop- 
erties of 7— matrices it is easy to show that only the Lowest Landau Level (LLL) 
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gives a non-vanishing contribution to the trace: 



L f =N c 



Qf\eB\ r dp 3 m q 
2% J 2% Ci)o(P3) 



(29) 



where 0)o = <£>k=Q- 

From Eq. d27b we notice that the prescription ( fT3l l is almost equivalent to the 
introduction of a running effective quark mass, 



that can be considered as a rough approximation to the effective running quark mass 
in QCD Wf\ which decays at large quark momenta, see also the discussion in |96l . 
Once the scale A is fixed, the Landau levels with n > 1 are removed from the chiral 
condensate if eB A 2 . 

In the upper panel of Fig. [5] we plot the chiral condensates for u and d quarks, 
as a function of eB, for the QM model. The magnetic field splits the two quantities 
because of the different charge for the two quarks. The small oscillations, which are 
more evident for the case of the M-quark, are an artifact of the regularization scheme, 
and disappear if smoother regulators are used, see the discussion in ||38l . In the 
regime of weak fields, our data are consistent with the scaling (//) <*; \eB\ 2 /M where 
M denotes some mass scale; in the strong field limit we find instead (//) ~ \eB\ 3 / 2 . 
The behavior of the quark condensate as a function of magnetic field is in agreement 
with the magnetic catalysis scenario. 

In the middle panel of Fig. [5] we plot our data for the expectation value of the 
magnetic moment. At weak fields, Ef °< \eB\ as expected from Eq. ([28). In the strong 
field limit, non-linearity arises because of the scaling of quark mass (or chiral con- 
densate); we find Ef |eB| 3 / 2 in this limit. 

In the lower panel of Fig. [5] we plot our results for the polarization. Data are 
obtained by the previous ones, using the definition d24l i. At small fields, the polar- 
ization clearly grows linearly with the magnetic field. This is a natural consequence 
of the linear behavior of the magnetic moment as a function of eB for small fields, 
see Fig. [3] On the other hand, within the chiral models we measure a saturation of 
Hf at large values of eB, to an asymptotic value jttoo = 1. This conclusion remains 
unchanged if we consider the NJL model, and it is in agreement with the recent 
Lattice findings l92l . It should be noticed that, at least for the u— quark, saturation 
is achieved before the expected threshold for p— meson condensation [98 99 100]. 
Therefore, our expectation is that our result is stable also if vector meson condensa- 
tion is considered. 

The saturation to the asymptotic value /Xo» = 1 of polarization is naturally under- 
stood within the models we investigate, as a LLL dominance in the chiral condensate 
(i.e., full polarization). As a matter of fact, Ef and (//) turn out to be proportional 
in the strong field limit, since only the LLL gives a contribution to the the latter, 
comparing Eq. ( f27T > and d29l which imply 




(30) 
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Fig. 5 Upper panel. Chiral condensates of u-quarks (red) and rf-quarks (blue), in units of the 
same quantities at zero magnetic field, as a function of the magnetic field. Middle panel. Expec- 
tation value of the magnetic moment operator, in units of /J , as a function of eB. Lower panel. 
Polarization of u-quarks (red) and d— quarks (blue) as a function the magnetic field strength. From 
Ref. 1951 . Copyright(2012) by American Physical Society. 



, (//)hll ,_ n 

WW ' (31) 

where (//)hll corresponds to the higher Landau levels contribution to the chiral 
condensate. In the strong field limit (//)hll — > 0, see Eq. (127k hence, jj.f has to 
approach the asymptotic value ji^ = I. On the other hand, in the weak field limit 
(//)hll (//) and the proportionality among £ f and (//) is lost. 
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At small fields fif = \%QfeB\ from Eq. d22l) . Hence, we use the data on polariza- 
tion at small fields, to obtain the numerical value of the magnetic susceptibility of 
the chiral condensate. Our results are as follows: 

% « -4.3 GeV 2 , NJL (32) 
X ~ -5.25 GeV~ 2 , QM (33) 

respectively for the NJL model and the QM model. To obtain the numerical values 
above we have used data for eB up to 5m 2 ps 0.1 GeV 2 , which are then fit using a 
linear law. Using the numerical values of the chiral condensate in the two models, 
we obtain 

*<//)« 69 MeV, NJL (34) 
*<//)« 65 MeV, QM (35) 

The numerical values of X that we obtain within the effective models are in fair 
agreement with recent results, see Table I. In our model calculations, the role of 
the renormalization scale is played approximately by the ultraviolet cutoff, which is 
equal to 0.560 GeV in the QM model, and 0.627 GeV in the NJL model. 

To facilitate the comparison with previous estimates, we review briefly the frame- 
works in which the results in Table I are obtained. In l87l the following result is 
found, within OPE combined with Pion Dominance (we follow when possible the 
notation used in |92l ): 

pd N c 

X = — Cv - 5 - , Pion Dominance (36) 

with F n = V2f n = 130.7 MeV and c x = 2; the estimate of ll87l is done at a renor- 
malization point M = 0.5 GeV. It is remarkable that Eq. d36b has been reproduced 
recently within AdS/QCD approach in [89|. Probably, this is the result more compa- 
rable to our estimate, because the reference scales in l87l and in this article are very 
close. Within our model calculations we find c% JL = 1.93 and cf 4 = 2.36. Using 
the numerical value of F n and c x we get x PD = —4.45 GeV -2 , which agrees within 
3% with our NJL model result, and within 18% with our QM model result. 

In 1881 the authors find c x =2.15 within hard- wall holographic approach, at the 
scale M< 1 GeV. The results of [88] are thus in very good parametric agreement 
with l87l : on the other hand, the numerical value of F K in the holographic model 
is smaller than the one used in [87], pushing the holographic prediction for x to 
slightly higher values than in f87l . However, the scale at which the result of 1 88 1 is 
valid should be much smaller than M = 1 GeV, thus some quantitative disagreement 
with f 87j is expected. As the authors have explained, it might be possible to tune the 
parameters of the holographic model, mainly the chiral condensate, to reproduce 
the correct value of F K ; their numerical tests suggest that by changing the ratio 
(/ f) l m p of a factor of 8, then the numerical value of c x is influenced only by a 5%. 
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Method 


X (GeV--) 


Ren. Point (GeV) 


Ref. 


Sum rules 


-8.6±0.24 


1 


|82| 


Sum rules 


-5.7 


0.5 


1 83 1 


Sum rules 


-4.4±0.4 


1 


|84| 


Sum rules 


-3.15±0.3 


1 


|85| 


Sum rules 


-2.85 ±0.5 


1 


I86| 


OPE + Pion Dominance 


-N C /(4^F£) 


0.5 


|87| 


Holography 


-1.075W c /(4^Fi) 


< 1 


1 88 1 


Holography 


-N c /(4n'F£) 


< 1 


I89| 


Instanton vacuum 


-2.5±0.15 


1 


|90| 


Zero mode of Dirac Operator 


-3.52 


1 


1911 


Lattice 


-1.547(3) 


2 


I92| 


NJL model 


-4.3 


0.63 


This work 


QM model 


-5.25 


0.56 


This work 



Table 2 Magnetic susceptibility of the quark condensate obtained within several theoretical ap- 
proaches. In the table, F K = 130.7 MeV. See the text for more details. Adapted from Ref. |95| . 
Copyright(20 1 2) by American Physical Society. 



It is therefore plausible that a best tuning makes the quantitative prediction of |[88l 
much closer to the estimate of 11871 . 

In l90l an estimate of % within the instanton vacuum model has been performed 
beyond the chiral limit, both for light and for strange quarks (the result quoted in 
Table I corresponds to the light quarks; for the strange quark, Xs/Xu,d ~ 0.15 is 
found). Taking into account the numerical value of the chiral condensate in the in- 
stanton vacuum, the numerical estimate of l90l leads to % = — 2. 5 ±0.15 GeV -2 
at the scale M = 1 GeV. An analytic estimate within a similar framework has been 
obtained in 11911 . in which the zero-mode of the Dirac operator in the background 
of a SU (2) instanton is used to compute the relevant expectation values. The result 
of US gives % = -3.52 GeV 2 atM « 1 GeV. 

In [92 1 the result % = —1.547 GeV~ 2 is achieved within a two-color simulation 
with quenched fermions. It is interesting that in [92] the same quantity has been 
computed also at finite temperature in the confinement phase, at T = 0.82T C , and 
the result seems to be independent on temperature. The reference scale of ||92l , 
determined by the inverse lattice spacing, is M sa 2 GeV. Therefore the lattice results 
are not quantitatively comparable with our model calculation. However, they share 
an important feature with the results presented here, namely the saturation of the 
polarization at large values of the magnetic field. Finally, estimates of the magnetic 
susceptibility of the chiral condensate by means of several QCD sum rules exist [ 82 , 
[83] El [85] [Ml. The results are collected in Table I. 
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In this Section, we make semi-analytic estimates of the polarization and the mag- 
netic susceptibility of the quark condensate, as well as for the chiral condensate in 
magnetic background, within the renormalized QM model. This is done with the 
scope to compare the predictions of the renormalized model with those of the effec- 
tive models, in which an ultraviolet cutoff is introduced to mimic the QCD effective 
quark mass. 

In the renormalized model, we allow the effective quark mass to be a constant in 
the whole range of momenta, which is different from what happens in QCD |[97ll . 
Thus, the higher Landau levels give a finite contribution to the vacuum chiral con- 
densate even at very strong fields. This is easy to understand: the ultraviolet cutoff, 
A, in the renormalized model can be taken larger than any other mass scale, in par- 
ticular A > \eB\ l l 2 ; as a consequence, the condition p 2 + 2n\eB\ < A 2 is satisfied 
taking into account many Landau levels even at very large eB. The contribution of 
the higher Landau levels, once renormalized, appears in the physical quantities to 
which we are interested here, in particular in the chiral condensate. 

Since the computation is a little bit lengthy, it is useful to anticipate its several 
steps: firstly we perform regularization, and then renormalization, of the QEP at 
zero magnetic field (the corrections due to the magnetic field turn out to be free of 
ultraviolet divergences). Secondly, we solve analytically the gap equation for the a 
condensate in the limit of weak fields, and semi-analytically in the opposite limit. 
The field-induced corrections to the QEP and to the solution of the gap equation are 
divergence-free in agreement with 1321 . and are therefore independent on the renor- 
malization scheme adopted. Then, we compute the renormalized and self-consistent 
values of the chiral condensate and of the magnetic moment, as a function of eB, 
using the results for the gap equation. Within this theoretical framework, it is much 
more convenient to compute (//) and Ef by taking derivatives of the renormalized 
potential; in fact, the computation of the traces of the propagator in the renormalized 
model is much more involved if compared to the situation of the non-renormalized 
models, since in the former a non-perturbative (and non-trivial) renormalization pro- 
cedure of composite local operators is required [?]. Finally, we estimate %, as well 
as the behavior of the polarization as a function of eB. 

5.2.1 Renormalization of the QEP 

To begin with, we need to regularize the one-loop fermion contribution, namely 

VfS£ = -NcZ^LPn r^(k 2 + 2n\Q f eB\+ m 2 q f 2 . (37) 

/ L% n=0 J -°° L% 

To this end, we define the function, y{s), of a complex variable, s, as 
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m=-Ncl}^t^r^ + 2n\Q f eB\+rn^ . (38) 



2k „ J-oo 2k 

The function ^{s) can be analytically continued to s = 0. We define then Vf"™™ = 
lim s ^ + 'f(s). After elementary integration over k, summation over n and taking the 
limit s — >• + , we obtain the result 

« P " =NcZ^^ -log(2|e / «B|) +«) B 2 (g) 

-««E^0-io S K) + «), 

where we have subtracted terms which do not depend explicitly on the conden- 
sate. In the above equation, is the Hurwitz zeta function; for Re(f) > 1 
and Re(q) > 0, it is defined by the series C(?,g) = LJT=o( w + <?)~'> m e series 
can be analytically continued to a meromorphic function defined in the complex 
plane t ^ 1. Moreover we have defined q = (ra;J + 2 1 Q/^-B | ) /2 1 g/-e.B | ; furthermore, 
a = 1 — Ye — y/(— 1/2), where y# is the Eulero-Mascheroni number and \j/ is the 
digamma function. The derivative (— l,q) = d£(t,q)/dt is understood to be com- 
puted at t = — 1 . 

The first two addenda in Eq. d39t arise from the higher Landau levels; on the 
other hand, the last addendum is the contribution of the LLL. The function B2 is 
the second Bernoulli polynomial; using its explicit form, it is easy to show that the 
divergence in the LLL term in Eq. d39l l is canceled by the analogous divergence in 
the first addendum of the same equation. It is interesting that the LLL contribution, 
which is in principle divergent, combines with a part of the contribution of the higher 
Landau levels, leading to a finite result. This can be interpreted as a renormalization 
of the LLL contribution. On the other hand, the remaining part arising from the 
higher Landau levels is still divergent; this divergence survives in the B —> limit, 
and is due to the usual divergence of the vacuum contribution. We then have 

v£gg=NcZ^i Q-log(2|<2 /e £|)+a 

/ 

Ar ^ \Q f eB\m 2 m 2 
+N ^^K 2 ~^m^B\ 

-*cL^C(-l,9). (40) 
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The renormalization procedure of the quantum effective potential is discussed in 
some detail in |95l . Here it is not necessary to discuss this procedure, and we just 
focus on the results. 



5.2.2 Approximate solutions of the gap equation 

Weak fields. In the weak field limit (eB <C m 2 ) the correction due to the magnetic 
field to the quantum effective potential can be computed: 

V, « -N C V log = —N c Y log < (41) 

1 Ll f 24tt 2 g 2\Q f eB\ cl f 24k 2 g jU 2 ' 

which is in agreement with the result of 1 32 ] . In the above equation we have followed 
the notation of |94| introducing an infrared scale ji, isolating and then subtracting 
the term which does not depend on the condensate. The scale fi is arbitrary, and we 
cannot determine it from first principles; on the other hand, it is irrelevant for the 
determination of the a— condensate. We expect H~f n since this is the typical scale 
of chiral symmetry breaking in the model for the a field. 

In this limit, it is easy to obtain analytically the behavior of the constituent quark 
mass as a function of eB. As a matter of fact, we can expand the derivative of the 
QEP with respect to a, around the solution at B = 0, writing (a) = f K + 8a. Then, 
a straightforward evaluation leads to 

As anticipated, the scale jj. is absent in the solution of the gap equation. 

Strong fields. In the limit eB 3> m 2 , we can find an asymptotic representation of V\ 
by using the expansion £'(— l,q) = cq +c\(q— 1) valid for with cq = —0.17 

and c\ = —0.42. Then we find 

„ ml ( ml , \ 2\Q f eB\ ^ \Q f eB\ml 
* - -*cl^ [-i + \QfeB\ J log^-^I^^ct , (43) 

where we have subtracted condensate-independent terms. 

In the strong field limit it is not easy to find analytically an asymptotic repre- 
sentation for the sigma condensate as a function of eB; therefore we solve the gap 
equation numerically, and then fit data with a convenient analytic form as follows: 

mq = b\eB\ l ' 2 + C Jf (44) 
\eB\ 

where b = 0.32 and c = 32.78. At large fields the quark mass grows as |eB| 1//2 as 
expected by dimensional analysis; this is a check of the equations that we use. 
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5.2.3 Evaluation of chiral condensate and magnetic moment 

Chiral condensate. To compute the chiral condensate we follow a standard proce- 
dure: we introduce source term for //, namely a bare quark mass mr, then take 
derivative of the effective potential with respect to nif evaluated at mt — 0. For the 
weak field case we obtain 



1Z71 A m q 



On the other hand, in the strong field limit we have 

, - , N c m a ,, , 2\Q f eB\ 

Using Equations d42l i and ( 144) . we show that the chiral condensate scales as a + 
b(eB) 2 for small fields, and as |eB| 3 / 2 for large fields. 

Magnetic moment. Next we turn to the computation of the expectation value of 
the magnetic moment. The expression in terms of Landau levels is given by Eq. (29} , 
which clearly shows that this quantity has a log-type divergence. In order to avoid a 
complicated renormalization procedure of a local composite operator, we notice that 
it is enough to take the minus derivative of V\ with respect to B to get magnetization, 
^# l94| , then multiply by 2m/ Qf to get the magnetic moment. This procedure is 
very cheap, since the fi— dependent contributions to the effective potential are finite, 
and the resulting expectation value will turn out to be finite as well (that is, already 
renormalized). 

In the case of weak fields, from Eq. d4Tb we find 

Q f \eB\m a m 2 

On the other hand, in the strong field limit we get from Eq. (03]) 

ml 2\Q f eB\ 
Z f =N c -^log-^ r ^ . (48) 
47T Z mi 



The above result is in parametric agreement with the estimate of magnetization 
in l94l . In fact, m 2 m \eB\ in the strong field limit, which leads to a magnetization 
^C^BlogB. 

Using the expansions for the sigma condensate at small and large values of the 
magnetic field strength, we argue that Ef « \eB\ in a weak field, and Ef ~ |eB| 3 / 2 
in a strong field. 
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5.2.4 Computation of chiral magnetization and polarization 

We can now estimate the magnetic susceptibility of the quark condensate and the 
polarization as a function of eB. For the former, we need to know the behavior of 
the magnetic moment for weak fields. From Eq. d47b and from the definition ( f22T > 
we read 

,- , N c m a m 2 
Xiff) = "Tzr^logTrf =f(P) ■ (49) 



The presence of the infrared scale ji makes the numerical estimate of % uncertain; 
however, taking for it a value jj, « f K , which is the typical scale of chiral symmetry 
breaking, we have % (ff) ~ 44 Me V, which is in agreement with the expected value, 
see Eq. d23l ). 

Next we turn to the polarization. For weak fields we find trivially a linear depen- 
dence of jif on \QteB\, with slope given by the absolute value of % in Eq. j49l . On 
the other hand, in the strong field limit we find, according to Eq. (08]), 



m 



'i 



m i _ 12/1 

\Q f eB\ ~ 1 ~ b+\Q f \ 



M/^ -l 1 i7i „ D i ^ 1 ~ l~TF7~\ ' ( 50 ) 



where we have used Eq. (144) . This result shows that the polarization saturates at 
large values of eB, but the asymptotic value depends on the flavor charge. 

It is interesting to compare the result of the renormalized model with that of 
the effective models considered in the previous Section. In the former, the asymp- 
totic value of jif is flavor-dependent; in the latter, jif — > 1 independently on the 
value of the electric charge. Our interpretation of this difference is as follows: com- 
paring Eq. ( l50t with the general model expectation, Eq. PTI ). we recognize in the 
factor \Qf\/ (b+\Qf\) the contribution of the higher Landau levels at zero temper- 
ature, which turns out to be finite and non-zero after the renormalization proce- 
dure. This contribution is then transmitted to the physical quantities that we have 
computed. The trace of the higher Landau levels is implicit in the solution of the 
gap equation in the strong field limit, namely the factor b in Eq. (144) . and explicit 
in the additional \Qf \ dependence in Eq. (150) . A posteriori, this conclusion seems 
quite natural, because in the renormalization procedure we assume that the effec- 
tive quark mass is independent of quark momentum, thus there is no cut of the 
large momenta in the gap equation (and in the equation for polarization as well). 
In the effective models considered in the first part of this article, on the other hand, 
the cutoff procedure is equivalent to have a momentum-dependent effective quark 
mass, m q = gO& (A 2 — p 2 — 2n\Q feB | ), which naturally cuts off higher Landau lev- 
els when eB ^> A 2 . At the end of the days, the expulsion of the higher Landau levels 
from the chiral condensate makes /Xf — > 1 in the strong field limit. Our expectation 
is that if we allow the quark mass to run with momentum and decay rapidly at large 
momenta, mimicking the effective quark mass of QCD, higher Landau levels would 
be suppressed in the strong field limit, and the result d50b would tend to the result in 
Fig.H 
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6 Conclusions 

In this Chapter we have summarized our results for the phase structure of quark 
matter in a strong magnetic background. Our theoretical investigative tools are chi- 
ral quark models improved with the Polyakov loop, which allow to study simulta- 
neously chiral symmetry breaking and deconfinement. 

The main motivation of this series of studies is of a phenomenological nature. 
In fact, it has been shown that huge magnetic fields are produced in non-central 
heavy ion collisions. These fields might trigger the P— and CP— odd process dubbed 
Chiral Magnetic Effect (CME). Therefore, in order to make quantitative estimates 
of the observables which are sensitive to the CME, it is extremely important to 
understand how hot quark matter behaves under the influence of a strong magnetic 
background. The other side of the coin is that simulations show these huge fields 
have a very short lifetime: therefore the present studies should take into account this 
time dependence. Moreover electric fields, which we have neglected so far, are also 
produced in the collisions. Finally, the electromagnetic fields considerably depend 
on space coordinates on the scale of the volume of the expanding fireball. This 
dependence has been ignored in our studies, since the magnetic background is taken 
to be homogeneous in space and constant in time. 

Our results support the scenario of magnetic catalysis, which manifests itself in 
both an increase of the chiral condensate at zero temperature, and an increase of 
the critical temperature for chiral symmetry restoration. Moreover, depending on 
the interaction used, deconfinement may occur either together with chiral symme- 
try restoration, or anticipate it. The latter possibility, even if more fascinating than 
the former since it opens a window for the Constituent Quark Phase, seems to be 
excluded by lattice simulations. 

Recent lattice simulations show that the critical temperature for chiral symme- 
try restoration, T c , is strongly affected by the quark mass. In particular, for small 
quark masses (hence, for the u and d quarks) the critical temperature decreases with 
the magnetic field strength; on the other hand, T c increases with the magnetic field 
strength for the s quark. As we have discussed in the main body of this Chapter, 
it seems that self-consistent computations within chiral quark models are not able 
to reproduce this feature, even when quantum fluctuations are taken into account. 
Thus, it remains an open problem to understand this unexpected behavior of T c . Cer- 
tainly independent simulations performed by other groups are necessary to confirm 
the present results. 

We have also briefly summarized a computation of the magnetic susceptibility of 
the chiral condensate, %, and of quark polarization, fXf at zero temperature, based on 
the quark-meson model. The computed value of % is m agreement with most of the 
previous estimates, and with experimental data. Moreover, this model gives a simple 
interpretation of the saturation of jif observed on the lattice: at very large magnetic 
field strength, the quarks occupy the lowest Landau level, expelling higher levels 
from the chiral condensate; hence, chiral condensate turns out to be proportional to 
the quark magnetic moment, making the ratio (that is, polarization) just a constant. 
In the case of the non-renormalized model, this constant turns out to be equal to one 
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and flavor independent; on the other hand, in the case of the renormalized model, 
the constant is flavor dependent. The latter result is easily understood: the renor- 
malization procedure of the momentum indepentent interaction of the quark-meson 
model brings all the Landau levels into the renormalized chiral condensate. We ex- 
pect that the replacement of the simple interaction discussed here with a non-local 
one, which should mimick the quark self-energy measured on the lattice, will make 
the expulsion of the higher Landau levels active also in the renormalized model, 
hence reproducing the results of the non-renormalized model. 
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